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Abstract 

Drift-diffusion computer simulations model available in Synopsys' Sen- 
taurus TCAD is used to study electrical, I-V, and optical, I-L, character- 
istics of separate-confinement heterostructure laser based on AlGaAs. We 
investigate the role of the width and depth of Quantum Well (QW) ac- 
tive region, below and above the lasing threshold. The device properties 
depend on both, the number of bound QW states and on closeness of the 
highest bound states to conduction or valence band offset. The lasing 
action may not exist at certain widths or hights of QW, and the thresh- 
old current is a discontinuous function of these parameters. The effects 
are more pronounced at low temperatures. Discontinuities in characteris- 
tics are found, at certain conditions, in temperature dependencies as well. 
The carriers scattering time on QW is shown to have the crucial role on 
amplitude of discontinuities. 

1 Introduction 

Computer modelling of electronic devices is a relatively new approach towards 
study of physical phenomena occurring there as well optimizing their technical 
characteristics. Methodologicaly, this field of scientific and engineering activity 
may be placed between theory and experiment, not belonging however to ei- 
ther of them: persuing research of that kind requires theoretical understanding 
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of physics of microscopic processes and may work as a helpfull tool in intcrp- 
tretation of experimental data. Contrary to sometime met thinking, computer 
modelling can not replace theory or experiment. In some situations, results of 
that research may provide an inspiration for understanding or testing physical 
phenomena: it is easier, faster and less expensive to perform modelling than 
experiments, and we are not restricted that much by, often large, inaccuracy of 
experimental data that may hide insightfuU details. 

When performing modelling AlGaAs SCH lasers with Synopsys' Sentaurus 
TCAD , we noticed unexpected steps in some of their characteristics (thresh- 
old current Ith versus the width of quantum well da (0), or versus it's height, 
etc). Analyses of results led us to conclusion that observed discontinuities occur 
when the most upper bound QW state crosses the conduction or valance band 
offset energy. Following that idea, we guessed that the effects may manifest 
itself in temperature dependence of other physical quantities as well, if laser 
parameters are choosen for that properly. In this work we show that the dis- 
continuities are found also below the lasing threshold current, in their I — V 
characteristics or in gain (or loss) versus current. This article is continuation of 
our efford to understand better the nature of these physical phenomena ([3 ). 
Here we show results of computer modelling of the role of quantum well (QW) 
scattering times. 

2 Modelling 

The laser we are modelling has dimensions, structure and doping as described by 
Andreev, et al. [1], [S]. The lasing wavelength is SOStito, the lasing offset voltage 
Uq is 1.56 — 1.60y, differential resistance just above the lasing offset, r ~ dU/dl, 
is 50 — SOmil, threshold current Ith is 200 — 300m^, slope of optical power, 
S = dL/dl, is 1.15 — 1.25W/A, and left and right mirror reflection coefficients Ri 
and Rr are 0.05 and 0.95. The reference laser has the width of QW, da, of 12nm 
and both waveguides' width is 0.2/im. In order to reproduce laser characteristics 
in computational results, we played with several variables available in Synopsys. 
The critical one is Aph - the effective surface area factor in Physics section of 
Synopsys command file. An agreement between experiment and calculation is 
reached for Aph of about 0.059. This low value of Aph should not be surprising: 
the physics we are dealing with here, in particular in waveguide and in QW 
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regions, depends on ballistic transport as well, while the computational model we 
use is derived from drift-diffusion equations, modified for dealing with transport 
through QW as discussed in the next section. 
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Table 1: A few sets of simulation conditions (A — F) for data shown in Figures 
andlH C is radiative recombination rate (Eq. [2]). a„, and ap are coefficients of 
free carrier absorption formula ([1]). Temperature for all cases is T = 300K, elec- 
tron and hole scattering times are assumed 8- 10"i'^s and 4- IQ-i'^s, respectively, 
and electron and hole mobility 9200cm^ /Vs and iOOcm^ /Vs, respectively. No 
additional light scattering mechanisms are considered. 

Other parameters available in Synopsys, important in this case, are these 
related to light absorption and carrier scattering. Experiment shows that ab- 
sorption coefficient is of the order of lcm~^ ([4')- It is argued that in AlGaAs 
lasers the main contribution to light absorption is due to photon scattering on 
free carriers, with the free carrier absorption coefficient, ctfc, given by: 

afc = {ctn-n + ap-p) ■ L, (1) 

where n and p are the electron and hole density, and L is light intensity. We 
choose in our calculations such values of a„ and ap that an effective absorption 
coefficient obtained would be close to that experimental one. 

It is important also to have a reasonable value of radiative recombination 
rate, Rr^ which is assumed to be described by: 

Rr=C-(n-p-nl^^), (2) 

where rii^^^ describe the effective intrinsic density, and C is a parameter 
available for changes. 

Typical I — V characteristics computed at T = 300 are shown in Figure [I] 
for a broad range of QW widths. For current near the lasing threshold current 
Ith (i-e. for voltage near the lasing offset voltage Ua), which correspond to a 
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Figure 1: Typical I — V characteristics computed at T = 300. OpticalLoss 
parameter is assumed 0, no radiative recombination, free carrier scattering rate 
parameters t„ and Tp are 8 • 10~^'^s~^and 4 • 10~'^'^s~'^, with electron and hole 
mobolities 9200cm^ /Vs and AOOcm^ /Vs. The legend describes width of QW 
(in nm), from 5nm from right-bottom curve to 2Anm for uppermost curve. 

kink in I — V, for most of these curves the results are very well approximated 
by a phcnomcnological modified exponential relation ([6]): 

IiU) = Ith-exp{A-{U-Uo)+B-{U^Uo)^), forU<Uo 
IiU)^Ith-exp{C-{U~Uo)+D-{U-Uo)^), forU>Uo ^' 

where Ith, Uq, as well A, B, C, and D are certain fiting parameters. 

There straightforward interpretation of these curves around the lasing thresh- 
old, since a strong, nonlinear interplay between the effects of carrier transport 
and scattering takes place, with light absorption as well. However, we may no- 
tice an interesting feature for parts of curves below the lasing threshold. While 
the width of QW, da, changes (nearly) monotonically, the curves however are 
grouped into a few sets such that they nearly coincide together, within each 
group. 

A very similar feature is observed when gain or loss is drawn as a function 
of current, for many widths of active region (Figure Again, we see that 
datacurves for current below the lasing threshold (that corresponds to kinks 
in curves) are grouped into a few sets such that they nearly coincide together 
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Figure 2: Current as a function of QW width derived at constant voltage from 
datacurves as these shown in Figure [1] and described in Table [T] The solid line 
for datapoints F is to guide the eyes, only. F is computed for constant Al 
concentration in QW of 8%, while all other datasets {A — E) are computed with 
such a concentration of Al in QW that lasing wavelength will remain constant 
(SOSnm) when QW width changes. The arrows are for datapoints F. From 
results of a separate analysis, not shown in this Figure, it follows that the 
relative hight of steps and position of these steps does not depend on voltage at 
which current is messured, even though current values may change as much as 
100 times. 
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Figure 4: Lasing threshold current as a function of QW width for datasets as 
these shown in Figure [T] and described in Table [TJ The arrows are at positions 
close to but not identical to these in Figured] 



within each group. If gain or loss were drawn as a function of voltage, however, 
we would not see such a grouping. 

Therefore, we conclude that below the lasing threshold, current as a function 
of QW width at constant voltage derived from data like these in Fig. [TJ or gain 
or loss as a function of QW width, at constant current values, also below the 
lasing threshold, will follow step-like functions. 

This is illustrated in Fig. [21 where current as a function of QW width derived 
at constant voltage from datacurves similar to these as in Figure [1] is shown. 
Figure [5| presents data computed at different conditions, and marked from A 
to F, for several combinations of free carrier scattering coefficients, q;„ and Op, 
and values of C, the radiative recombination parameter, as described in Table [T] 
Moreover, dataset F differs from datasets A-E. The last are computed assuming 
changing Al concentration in waveguides (when Al in QW is kept constant) in 
such a way that the lasing wavelength does not change with the change of QW 
width. The dataset F is computed for constant Al concentration in waveguides 
of 33%. The solid line in Fig. [5] is drawn through datapoints F, and arrows 
there refer to curve F as well, and indicate positions of bound QW energy states 
crossing the conduction- or valance band offset energies. Positions of bound QW 
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Figure 5: T=300K. Conduction and valence band offset energy, {Ecbo and 
EvBo), and electron (En), light- and heavy hole energies (LHn and HHn) in 
QW, as a function of QW width. Hole energies and Evbo have been scaled up 
by a factor 28 to obtain coincidnce with electron energy scale (i.e., Ecbo and 
Evbo curves are the same in this Figure). 



energy states for cases A-E, as illustrated in Fig. [SJ are very close to but not 
identical. The results of Figure [5j described in more details earlier ([7]), were 
found also to be exactly the same as computed by using nextnano software ([8 ). 

The step- like features are preserved also in Ith{da) dependencies, as illus- 
trated in Fig. 01 There, however, effects of carrier scattering and light absorp- 
tion smear-out the picture. It is useful to notice that at some values of QW 
width no lasing action is reached, and therefore the datapoints in that Figure 
are not available for all QW widths. 

Changes of QW height (caused by difference of Al concentration in QW and 
waveguides) cause very similar step-like dependencies. Moreover, the effects are 
in some situations more clear and pronounced at low temperatures. We did 
modelling for T=77.6 K to confirm their existence. 

Moreover, with a careful design of laser structure (content of Al in QW 
and waveguides) it is possible to find the evidence of the effect in temperature 
dependence of current, when measurements are performed at constant voltage. 
We found such Al concentrations when the number of QW bound states changes 
with temperature swap. Figure [5] shows how the uppermost bound electron state 
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Figure 6: The energy difference between conduction band ofFset energy Ecbo 
and the closest to it bound electron state in QW, £^4, as a function of tempera- 
ture, when active region concentration is 8% of Al, for three values of waveguide 
Al concentraion: 34.60%, 34.65%, and 34.70%. 



energy, in this case, differs from Ecbo (energy levels of other QW bound 
states do not play a significant role in this case), for a three Al concentrations in 
waveguide, when Al concentration in QW is 8%. For Al concentration 34.70%, 
the -E4 energy level exists always through temperature swap studied. For Al 
concentration 34.65%, the £4 energy level does not exist between around 150 and 
370A'. For Al concentration 34.60%, it exists at temperatures higher than about 
470Ar, only. This has profound implications on I{T) dependencies measured at 
constant voltage for these three different Al concentrations of Al in waveguide, as 
Figure[7]illustrates. For 34.60% and 34.70% of Al content, we observe continuous 
I{T) curves. However, for 34.65% of Al content, at low temperatures I{T) 
results fall on curve that has been computed for 34.60% of Al, and at high 
temperatures they fall on the curve computed for 34.70% of Al. 

3 The role of QW scattering times. 

The modeling results reported above were computed by using default (in Synop- 
sys TCAD) values of parameters related to carriers scattering times on QW (i.e. 
time of living of carriers in bound QW states) . An approach used in Synopsys 
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Figure 7: Temperature dependence of current for three values of voltage applied 
(1.510V, 1.613V, and 1.687V; for each of these this is below the lasing threshold), 
and for free values of Al concentrations in waveguide, the same as in Figure IHl 

is to have constant vaues of these parameters (which however may be adjusted 
by user). 

Existing theoretical and other computational results indicate that carriers 
scattering times depend on QW width ([9], [10]). Scattering on longitudinal op- 
tical phonons, LO, is considered as the main mechanism. Birner 10,) computed 
lifetime of electron states for transitions between the initial state E2 and the 
final ground state Ei), i.e. for an intersubband transition, for different quantum 
weU widths at T = 0. 

The nextnano^ calculations of Birner (flOJ) are in good agreement with these 
of Ferreira and G. Bastard ([S]). For quantum well widths smaller than about 
5.4 nm, only the ground state is confined and E2 is unbound. For quantum well 
widths larger than about 18nm, the transition energy becomes smaller than the 
LO phonon energy and scattering through the emission of an LO phonon is not 
possible any more. 

In order to have a better view on what is the effect of scattering times on 
the amplitude of steps (and their existence) in opto-electrical characteristics of 
the device, we performed simulations for several sets of values of parameters 
available in Synospsys, as described in Tabled] 
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Table 2: A few sets of simulation conditions {a — h in the first row) for data 
shown in Figures [9] - ITT] QW width is 12 nm. 



Typical results are shown in Figure \8\ Shorter QW scattering times lead 
to strong increase of current and pronounced steps in current as a function 
of QW width. It was found that the value of parameters eQWMobility and 
hQWMobility has no a noticable effect to these results. 

For a completness of the results on modeling the role of QW scattering, we 
show also I-V characteristics (Figure IH]), lasing light intensity L-I (Figure [TU)) 
and optical efficiency curves fFigurc fTTj) for a few sets of values of QW scattering 
parameters. 

4 Discussion 

In case of tunneling energy barrier, transfer matrix approach is used to describe 
charge transport through it ([14]. |15)). The interband tunneling current is 
written as 

- IeZ: •^(^) • fiE) ■ T{E) ■ dE (4) 
where T{E) is energy-dependent tunneling rate, N{E) is the density of 
states, f{E) is the Fermi-Dirac distribution function, respectively, and Emin 
and E,nax are minimum and maximum carrier energies available. 

In Scntaurus, s simplified intuitive model is used to handle the physics of 
carrier scattering at the quantum well (Figure [T^]). The carrier populations are 
separated into bound and continuum states, and separate continuity equations 
are applied to both populations. The QW scattering model accounts for the net 
capture rate, that is, not all of the carriers will be scattered into the bound states 
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da [nm] 



Figure 8: Height of steps does not depend in a noticable manner on parame- 
ters eQWMobility and hQWMobility. The data in this figure are ah computed 
by using the same values of these parameters, 9200 and 400 err? jV s, respec- 
tively. The effect of QW scattering times is however pronounced. The data, 
for curves A to D (from top to bottom) were computed with the following sets 
of (QWeScatTime, QWhScatTime): (10-l^ 5 • 10""), (8 • 10-i'\ 4 • IQ-"), 
(4 • 10-12^ 2 • 10~i2)^ aj^j (10"", 5 ■ IQ-^^)^ where time is in seconds. 
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Figure 9: I-V characteristics for a laser with 12 nm QW width, with QWeScat- 
Time, QWhScatTime, eQWMobihty, and hQWMobihty parameters, for curves 
from a to h, as described in Table [5J The curves with the same values of 
QWeScatTime and QWhScatTime coincide together, regardles of of values of 
eQWMobihty and hQWMobihty. 




0.5 1 1.5 2 2.5 3 3.5 4 

I [A] 

Figure 10: L-I characteristics for a laser with 12 nm QW width, for cases a to 
h the same as these in Figure M and described in Table [2] 



13 




0.5 1 1.5 2 2.5 3 3.5 4 

I [A] 

Figure 11: Optical efficiency, L/{U ■ I) for a laser with 12 nm QW width, for 
cases a to h the same as these in Figure [3] and described in Table [H 

of the quantum well. The electron capture rate from the continuum (subscript 
b) to the bound (subscript b) states is: 

i? = J^JE, J^JEt, ■ N,iE,) ■ NbiEb) ■ SiEb, E,) ■ f,iE,){l ~ MEb)) (5) 

where Ec and Ei, is energy of lowest conduction band-, and bound QW elec- 
tron states, N{E) is the density-of-states, S{Et,, Ec) is the scattering probability, 
and f{E) is the Fermi-Dirac distribution. The reverse process gives the electron 
emission rate from the bound to continuum states: 

M = /J^ dEc dEb ■ Nc{Ec) ■ m{Eb) ■ S{Eb, E^) ■ h{Eb){l - MEc)) (6) 

The net capture rate is C = R ~ M , and for very deep quantum wells 
(keyword QWDeep must be used for that in Sentaurus) is known to be given 
by approximation: 

C = R-M^{l-exp{T]b~Vc))-— (7) 

T 

where rji, — {—q^t ~ Ec)/kBT and rjc = {—q^c — Ec)/kBT contain the quasi- 
Fermi level information and r is the capture time. The capture time represents 
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Figure 12: The model of carrier scattering at the quantum well used in Sentaurus 
(based on [1]). 

scattering processes attributed to carrier-carrier and carrier-LO phonon inter- 
actions involving bound quantum well states (of which, it is generally assumed, 
carrier-LO phonon is dominating in the case considered). The net capture rate 
C is added to the continuity equations as a recombination term. 

In a similar way scattering of holes is treated, with their own characteristic 
capture time. These parameters are specified in Sentaurus by the keywords 
QWeScatTime and QWhScatTime. Their default values, 8 • 10~^^s and 4 ■ 
10~^^s, respectively, correspond reasonably well to these based on theory f[TT]. 
[12|). while photoluminescence spectroscopy results give values of an order of 
3 — 20ps [12) . In most of our modelling, if not indicated otherwise, we use 
also default values of electron and hole mobility, represented in Sentaurus by 
parameters eQWMobility = 9200cm'^/Vs and hQWMobility = AOOcm^/Vs, 
since we do not observe a noticablc changes of I-V-L laser characteristics when 
these parameters change a few times in any direction. 

For shallow quantum wells, the energy transfer during scattering can only 
occur in a limited range. In the limit of elastic scattering, the net capture rate 
is then approximated by: 
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\F,/M F,/MJ r ^ > 

where Fm is the Complete Fermi-Dirac integral of the order of m. The 
shallow quantum well model is activated by the keyword QW Shallow. 

It should be pointed out that Equations [7] and HI for deep and shallow quan- 
tum wells, respectively, while provide a convenient, intuitive description of car- 
riers scattering and capture on QW, these are approximate only. In particular, 
there is no dependence of capture time on energy of unbound carriers there and 
no periodic oscillations as a function of QW size ([H], [12] )• GaAs/AlGaAs are 
considered to have deep quantum wells. However, as we have shown, the results 
of our modelling indicate on a strong role of bound QW states located very 
closely to the offset energy levels of quantum wells. For these reasons, we did 
not restrict our calculations to deep- or shallow- QW models but used instead 
the full model available in Sentaurus. 

Equations S] - [5] all depend on density of bound states in QW. We expect 
hence that current through the QW will be proportional to the density of all 
bound states in QW. In effective mass approximation, the two-dimensional den- 
sity of electron states within each QW subband n equals 

N„{E)^^, forE>E.r, (9) 

Hence, the current should be proportional to the number of bound states x 
carrier mass. The quantity computed this way (with a certain multiplication 
factor) is represented by large circles in Figured] Though it must not be exact 
(for instance, no difference in scattering rates for electrons and holes is accounted 
for), it fits reasonably well I (da) dependence. 

In our modelling, we assumed scattering times independent of QW width 
(modeling is performed however for a few sets of scattering times). In Birner 
and Ferreira results ([5], [TO]), scattering time steeply increases with decrease of 
QW width, below the width of around 6 nm, and slowly, monotonically decreases 
when it becomes larger that about 6 nm, with a value of around lO^^^s at 7 
nm. 

The default value of scattering time used by Synopsys is 8 • 10~^^s, for 
electrons. That suggests that results represented by data on curves B and C in 
Figure [5] should be closest to these expected experimentally. At the same time. 
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changing of scattering times with QW width should not diminish the existence 
of steps since these changes are monotonic. Also, the steps should become more 
pronounced at large QW widths, as well the current values should increase in 
a steeper way than would follow from Figures [5] or [51 It ought however to be 
remembered that in Birners example results, only elecrons scattering is taken 
into account and only transitions between the lowest bound electron states. 
The results reported in literature often predict a quasi-periodic dependence of 
scattering times on QW width ([11], [H], [13). 

5 Summary and Conclusions 

When performing modeling of laser characteristics as a function of the width of 
active region we noticed a non-monotonic, discontinuous dependence of I (da) 
(when measured at constant voltage applied). A careful analysis of the data 
led us to the hypothesis that discontinuities occur when the most upper QW, 
bound energy states are found very close to the conduction or valence band 
energy offsets. The effect, hence, is thought to be related to changes in density 
of states of carriers from one hand, and to fast changes in carrier transfer matrix 
through QW for QW bound states close to Ecbo or Ey bo ■ As such, it ought to 
be more pronounced at lower temperatures, as confirmed by results of modeling 
I {da) at liquid Nitrogen temperature ([!]). 

The effect is observed also when modelling current as a function of QW 
depth (Al concentration in waveguide). 

Therefore, we concluded that a similar effect will be present also in modelling 
/ as a function of temperature. In that case however a carefull design of laser 
properties is needed, in such a way that a transition of the most upper QW 
energy state will pass through an edge of quantum well when temperature is 
swapt. 

We expect also that performing measurements on laser devices under uniax- 
ial or hydrostatic pressure might provide experimental evidence on signficance 
of these effects in real devices. 

These observations are potentially important for proper designing of semi- 
conducting lasers (choice of Al concentrations, thickness of the active region, 
etc), and, potentially, might be useful for designing a kind of quantum level 
spectroscopy tool for testing lasers for technological applications. 
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Comparison of results with these obtained when balhstic transport is in- 
cluded is very desirable. 
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